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different situations and have applications in several other branches of physicspT], 
[28]. Although solitons are well known and much studied, we believe that a simple 
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■^ I Abstract. A classification of the time evolution of the two-soliton solutions of the 

Boussinesq equation is given, based on the number of extrema of the wave. For solitons 

fT^ , moving in the same directions, three different scenarions are found, while it is shown 

that only one of these scenarios exists in case of oppositely moving solitons. 
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1. Introduction 

K* ■ Interest in environmental flows pr|-[T8] both utilize and motivate theoretical work in 

[^^ ■ nonlinear wave phenomena. Among these are solitons [19]- [26], which appear in many 

G\' 
l> 

^' 

O ' classification scheme of the evolving wave shapes of two colliding water wave solitons 

may not be in vain. 

Weakly nonlinear waves in a shallow channel may be well described in terms of the 
Boussinesq equation. It allows propagation in both directionsljl The small parameters 
Jh ■ of the problem are the ratio of the wave amplitude to the water depth (e = a/h) and 

" " the ratio of the water depth to the wave length {6 = {h/l)"^), these are assumed to be 

comparable. The smaller these parameters are, the better the approximation will be. 
Therefore, any phenomenon found when solving the Boussinesq equation corresponds 
to observable effects if e and 6 are sufficiently small. 

It is well known that the Boussinesq equation is fully integrable and has soliton 
solutions. Especially, when two solitons have collided, their individual properties 
are preserved and completely restored when they spatially separate. The change of 
the wave's shape during collisions is interesting and, as it turns out, follows three 
possible scenarios. Especially, it is intriguing, what happens if solitons with very similar 
parameters collide, since it is well known, that the one-soliton solution is unique (up to 
a scale transformation), hence, two identical, albeit spatially separated sohtons do not 
form an exact solution. 

I Strictly speaking, the derivation of that equation does not allow counterpropagating waves. 
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In the present paper we classify the possible wave shape changes during the collision 
of two Boussinesq solitons. 

2. General properties of the two-soliton solution of the Boussinesq equation 

The Boussinesq equation in dimensionless form is given by 

Tltt -Vxx-3 (?7^) - T^xxx'x = . (1) 

\ / XX 

Two-soliton solutions are generated by the Zakharov-Shabat scheme ^9j, and may 
be written as 

d j h{l + q) + k2il + p) - ih + k2)a \ ,. 

"^ dx\ (l + g)(l+p)-a J ' ^ ^ 

where the new variables p and q are defined by 

_ 2kix-2u]it _ 2k2X-2uJ2t /n\ 

and the parameter a is given by 

a = 2 • (4j 

In these expressions ki and k2 are arbitrary positive parameters characterizing the 
individual solitons and 



u, = ±k,^l + Ak], (j = l,2) (5) 

with positive sign for a soliton propagating from left to right, and negative otherwise. 
Note that if ui ^ uj2, 

< a < 1 . (6) 

On the other hand, if Ui = uj2, we also have ki = k2 that corresponds to a one-soliton 
solution, not being of interest here. 

An important observation is that in the extreme cases when /ci, k2 ^ 1/2 or 
fci, k2^ 1/2 the parameter a depends only on the ratio 

k2 
Explicitly, we get for solitons moving in the same direction 

(l+,)2+'i(l_,)2 if fcl,fc2>| 

and for solitons moving oppositely 

[ if fci,A;2<i 

For the spatial derivative of the wave we have 

dx ((i + g)(i+p) -a) „=ij=o 



(10) 
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where the coefficients are given by 
Cio= -kl{a-lf 
cn= -kl{a-iy 
C20 = -k^a- 1) 

C21 = - {ki + k2) a{ki + /c2)^ - Ski + 3kik2 - Skj 
C22 = - kl{a - 1) 
C31 = {ki - k2) [a{2kl + Sfcifca + 2kl) - Skf - 3kik2 - 3kl 

C32 = (^2 - ki) \a{2kl + 5kik2 + 2kl) - Skj - 3kik2 - 3kl 

C41 = — ki 

C42 = - {ki + k2) a{ki + k2f' - ^kf + 3kik2 - 3ki. 

C43 = — k2 

C52 = ^1 

C53 = ^2 

All the other coefficients are zero. Again, in the extremes ki, ^2 ^ 1/2 or ^i? ^2 ^ 1/2 
a common factor kf can be pulled out of all the coefficients and the rest depends only 
on the ratio k2/ki. Accordingly, for both very small and very large wave numbers the 
behavior of the solitons is determined by this ratio, up to a scaling. 

For minima and maxima of the amplitude 7]{x, t) at a given time t we have 



n=l j=0 



. 



:ii) 



Additionally, due to the definitions ([5]) we have 

fci -fc2 _ J2{k2uii-kiLLj2)t ('12') 

Thus, extrema are given by the intersections of the curves f fTTj) and f lT2|) . The graph 
of Eq.(1T2l) is a power function with a time-dependent coefficient. In fact, time appears 
only here. As for the graph of Eq. dlip . we may prove some general properties on the 
basis of Eqs. lfTOl) . namely 



(i) Exchanging ki with k2 and Ui with U2 is equivalent with exchanging p with q. This 
follows directly from Eqs.([2])-(|4]). 

(ii) The curves can have at most three intersections with a line p = const, or q = const. 
Indeed, for a fixed p Eq. (TTTl) is a third order polynomial of q and vica versa. 

(iii) The transformation 

j9^(l-a)-, g->(l-a)- (13) 

p q 

leaves Eq. (fTT]) invariant. 

Indeed, direct substitution shows that under transformation (IT^ the expression 
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En=iE"=oCnjyg" ^ goes ovei into 






p3q3 



(14) 

n=l j=0 

In view of the definition of p and q, the symmetry (TT3l) means that the two-sohton 
solution is invariant with respect to a simultaneous spatial and temporal reflection, 
with respect to suitably chosen origins. Explicitly, transformation flT3|) is equivalent 
with 



X — )■ 2a;n — X , t — )■ 2tn — t 



(15) 



where 



Xo 



tn 



i (X'2 — Wi 

4 kiU2 - k2UJi 
1 k2-ki 



ln(l - a) 
ln(l - a) . 



(16) 



4 A;ia;2 - ^2^1 

The symmetry manifests itself in the log-log plots FigJH FigJHl FigJSl Fig|71 FigJH] 
as an inversion symmetry with respect to the symmetry point 

p = q = Vl — a . (17) 

(iv) Although p = 0, q = satisfies Eq.f fTTj) . in the first quadrant no curve starts from 
the origin. This is because for small p and q the first order terms dominate, and 
(for positive p and q) they are both negative. 

(v) For p — )■ 00 two asymptotes exist, namely, one at g = and another one at g = 1. 
Similarly, for g — )■ 00 we have asymptotes at p = and p = 1. This result can 
easily be obtained since e.g. for large values of p the terms containing the highest 
(third) power oi p dominate, i.e., 

k^p^q^ - klp^q = . (18) 

(vi) Near the above asymptotes the curve may be approximated by 



g 




p 



kf 



k^ 1 

y^(l — a)- if g — )■ 00 , p — )• 



- if 



(19) 

if p -> 00 , g -^ 1 (20) 

(21) 

g -)■ 00 , p -)■ 1 (22) 



This can be readily shown by taking into account the next-to-highest (second) 
power of the large variable. 

The right hand sides of Eqs. (IT9|) . fl2T|) are positive, thus the curve approaches the 
asymptote from above and from the right, respectively|§| The same is true for 



Hence they lie in the first quadrant of the (p, q) coordinate system. 
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Eqs. lpll) . (122]) in case of oppositely moving solitons {ujiUJ2 < 0). In contrast, for 
solitons moving in the same direction {(jJiUJ2 > 0) the right hand sides of Eqs.f l20p . 
(122|) can be both negative or positive, depending on the parameters ki, k2- 

The above properties allow us to explain the possible topologies of the curve f lTTl) . 




Figure 1. The zeros of Eq. (|TT|) for solitons moving in the same direction. Left paneh 
Type I. case (at fci = 1, k2 — 1.5), right paneh Type II. case (at fci = 1, fc2 = 1.7). 



Under transformation (TT3|) the point p = 1, q = oo goes over into p = 1 — a, q = 0, 
similarly, the point p = oo, q = 1 goes over into p = 0, q = 1 — a. Finally, the point 
p = 0, q = oo goes over into p = oo, q = 0. This allows only three possible topologies, 
according to the fact that the point p=l — a, q = should be continuously connected 
with either p = 0, q = l — a (FigHJleft panel), or p = oo, q = (FigJH right panel), or 
p = 1, q = oo. This last possibility is visualized again by the right panel of Fig{T]if p is 
exchanged with q, which, according to property 1 above, is equivalent with exchanging 
ki with /c2- Since these curves are intersections of a smooth surface with a plain, other 
possibilities are ruled out. We shall call the topology shown in the left panel of Fig|T] 
the Type I. case, and the topology shown in the right panel of FigJUthe Type II. case. 
The topology obtained from the Type II. case via exchanging the axes will be called the 
Type II. case. Note that the same transformation does not change the topology Type 
I. 

While all the three situations do occur for solitons moving in the same direction, 
the Type I. case never occurs for oppositely moving solitons. This can be shown by 
considering the intersections of the curve fITT]) with the p = q line, i.e., the zeros of 



where ^ = p 



= q and 
6 = 2- 



(fci + k2f 



h^ -\- h^ 

Polynomial fl23|) always has zeros at ,^ = and C, = ±"\/T 
last quadratic factor is 

a{ki + k2f 



(23) 

(24) 
The discriminant of the 



D = ¥-A{1- a) 



{kf + klf 



a — 



I2kik2{kl + kl) 
{ki + k2f 



(25) 
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It is a simple exercise to show that the last factor on the right hand side of Eg. (1251) is 
always negative in the case of oppositely moving solitons (cf. Eq.(|l])). Hence, in that 
case only a single positive root exists, while the case Type I. shown in the left panel of 
Fig{T] requires three positive roots. 

When changing the parameters fci, k2-, the curves with different topologies go over 
into each other. If one considers solitons moving into the same direction, it is sufficient 
to consider the case ki < ^2, since the parameters must not coincide and the opposite 
case (i.e. ki > k2) simply corresponds to the exchange of the axes. While the relative 
difference of parameters is sufficiently large, we have the topology Type II. If the relative 
difference is diminished, we get the topology Type I. The crossover between the two 
topologies if shown in left panel of Figj2j 



* 




1 1 




inl: 


\ / 




U^ 




"-^ 




/ vf- 






/ 




q 10-1- 




1 1 10' 




IK'^ 




/ ^\ ""'^ 




»t 




\^ W-i 








Itl-^ 





10'- 10"' JCP 10' IC^ 



I0-* w- 15'' 



Kf \tf 10= 



Figure 2. Crossover of the curves (ITT|) having different topologies. Left panel: 
crossover from Type I. to Type II. for solitons moving in the same direction (at fci = 1, 
^2 = 1.57912575), right panel: crossover from Type II. to Type II. for solitons moving 
in opposite directions {&i ki — 1, k2 ~ I) ■ 



One might wonder whether the middle branch of the curve touches both the other 
branches indeed at the same parameters. The answer is affirmative, and follows again 
from the symmetry (IT3l) . At the crossings of the branches one gets zero first derivatives 
in two independent directionq]]], hence, one has simultaneously 



f{p,q,ki,k2) 
df{p,q,ki,k2) 

dp 
df{p,q,ki,k2) 

dq 








(26) 
(27) 

(28) 



where 



f{p,q,ki,k2 



71=1 J =0 



(29) 



This is actually an intersection of a saddle with the tangent plane. 
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Applying now the transformation f ll3p to a crossing point, it is straightforward to show 
that at the transformed point Eqs. fl26l) - fl28l) are also satisfied. Indeed, the two crossings 
go over into each other under transformation ( TT3|i ^. 

In case of oppositely moving solitons, ki = k2 is certainly possible, and a transition 
from Type II. to Type II. occurs exactly when the two parameters coincide (cf. right 
panel of Figj2]). 



3. Time evolution of the extrema 



Extrema of the waves are the intersection points of Eqs. lTTT]) and (TT2!) . As time goes 
on, the coefficient of the power function ( IT2l) changes from zero to infinity, and the 
corresponding curvCJ "sweeps through" the curve flTTj) . According to the possible shapes 
of curve ( !TT|) we observe different scenarios, according to the number of maxima. As 
shown in FigsI3]and|ll in the Type I. case we always get three intersections, i.e., there 
are always two maxima and a minimum between them. 
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Figure 3. Type I. topology for intersections of tire curve (fTTj) and ([T2|) at parameters 
fci = 1.0 and ^2 = 1.5 for solitons moving in the same direction. Different colors of the 
straight lines correspond to different time instants. 



^ The transformation is obviously an involution, i.e., it is equal to its own inverse. 
+ On the log- log plots the graph of Eq. (IT^ is a straight line. 
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Figure 4. Type I. topology as seen in wave pattern for solitons moving in the same 
direction. Parameters and time instants identified by the colors are the same as in 
FigEl 



In contrast, in the Type II. and Type II. cases there are periods when only a single 
maximum exists. However, one has to distinguish here two subcases: 
a) If the middle branch is steep enough, namely, if 



(9 In g ^2 
(91np ki 



(30) 



at the symmetry point (1171) . then at certain times the two maxima reappear (Figsj5]and 
E]). This situation will be called Type II. a (or Type II. a). 



In q 




Figure 5. Type II. a topology for intersections of the curve (fTTj) and (|12[) at parameters 
fci — 1.0 and ^2 = 1.6 for solitons moving in the same direction. Different colors of the 
straight lines correspond to different time instants. 
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Figure 6. Type II. a topology as seen in wave pattern for solitons moving in the same 
direction. Parameters and time instants identified by the colors are the same as in 
FigEl 



In such a case the initially well separated solitons (two maxima) coalesc, the 
remnant of the smaller soliton being only a "drooping shoulder" at the front side. 
Further on, the "shoulder" moves towards the maximum, it rises and develops a second 
maximum. Thus a shallow valley is created on the top of the wave. At later times 
these events take place in reversed order: the rear bank of the valley goes down, the 
corresponding maximum disappears and becomes a drooping shoulder at the rear side, 
then it develops a maximum again behind the taller wave and the two solitons are again 
separated. 

b) If the steepness of middle branch is smaller than ^2/^1) no valley is created on 
the top of the wave, as shown in Figs. [3 and [HI We shall call this case Type II. b (or 
Type II. b). 



In q 




Figure 7. Type II. b topology for intersections of the curve ([TT|) and (IT^ at parameters 
fci — 1.0 and ^2 — 1.8 for solitons moving in the same direction. Different colors of the 
straight lines correspond to different time instants. 
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Figure 8. Type II. b topology as seen in wave pattern for solitons moving in the same 
direction. Parameters and time instants identified by the colors are the same as in 
FigEl 



Both Type II. a and Type II. b can be observed in case of solitons moving in the 
same direction. In contrast, for oppositely moving solitons, only the case Type II. b (or 
Type II. b) can exist. This can be proven on the basis of Eqs. (!30|) . ( !T7|) . Indeed, we have 
at the symmetry point 

Aklil - a) - a{k2 - hf + ^/T^^{Akl - a{ki + ^2)^) 



d\nq 



d\np 4A;|(1 - a) + a{k2 - k^)^ + ^r^^(4A;i - a{ki + fca)^) 

which is smaller than ^2/^1 for oppositely moving solitons, if k2 > ki. 



(31) 
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Figure 9. Type II. b topology for intersections of the curve (ITT|) and ([T^ at parameters 
fci = 1.0 and fc2 ~ 1-5 for solitons moving in opposite directions. Different colors of 
the straight lines correspond to different time instants. 
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Figure 10. Type II. b topology as seen in wave pattern for solitons moving in opposite 
directions. Parameters and time instants identified by the colors are the same as in 
FigEl 



4. Parameter space 

The scenarios described above are summarized in parameter space in FigsJTT] and [T2l 




Figure 11. Scenarios shown in parameter space for solitons moving in the same 
direction. Legend: red, dark red: Type I. (see Figl3]and FiglU, white: Type II. a (see 
Figl5]and Fig|6]), gray: Type II. a, green: Type II. b (see Figl7]and FigIS]), dark green: 
Type IT.b. 



For solitons moving in the same direction, at the border between Type II. a and 
Type II. b in parameter space (see Fig. [TH) the expression (13T1) is equal to k2/ki. 
This condition defines the border. As noted after Eg. (1101) . for both ki , A;2 -C 1 and 
ki , ^2 ^ 1 the expression ( 1311) depends only on the ratio ^2/^1, hence the border looks 
linear. In fact, its slope slightly differs for small and large k values. 

As discussed above, the border between between Type I. and Type II. a in parameter 
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space is given by Eqs.( EB]) -( EU]) . Again, the border is not exactly a straight hne, its slope 
is slightly different for large and small k values. 

The shaded and the unshaded regions in FigJTTl are obtained by exchanging ki and 
k^. 

As noted before, the fci = ^2 line is not allowed. 

For solitons moving in opposite directions, the parameter space is even simpler (see 
FigJT2l). All parameter values are allowed, and the crossover from Type II. b to Type 
II. b occurs at ki = k2. 




4 6 



Figure 12. Scenarios sliown in parameter space for solitons moving in opposite 
directions. Legend: green: Type II. b (see Figl9]) and Fig ITU] ), dark green: Type 
Il.b. 



5. Nearly identical solitons 



Let us consider now the situation when two nearly identical solitons interact. If they 
were strictly identical, we would obtain a single soliton solution rather than a two-soliton 
solution. Hence approaching the limit can be interesting. This situation corresponds to 
case I. Maxima correspond to intersections of the line (I12p with the two outer segments 
of the graph of (TTT]) . Obviously, as time goes on, the initially well separated solitons 
approach each other, then, without coalescing, their distance grows again. The minimal 
distance between them may be estimated (cf. Fig. flT3|) ) as 

({ki + k:^^ 



6t ■ > 



ki + k 



ln(l- 



ki + k 



In 



^ 12 Vfci 



i±!2 V 
k2j 



\{h 



ko 



J_ 

12 



fci 



U!2 
fc2 



.(32) 



Evidently, 5xmin diverges logarithmically as ki —> k2. 

If one considers now the corresponding waves (Fig. flT4l) ). initially one sees two very 
similar solitons, the (slightly) taller one chasing the smaller one. When approaching, 
the tall soliton gradually loses its height and speed, at the same time, the smaller soliton 
gains hight and speed, the distance between them starts increasing, and eventually, we 
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see the original solitons again, but this time the smaller one chasing and the taller one 
escaping. 
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Figure 13. Intersections of the curve (fTTj) and (fT2|) at parameters fci = 1.0 and 
^2 = 1-1 for solitons moving in the same direction. 






Ui. 



Ill 



LA. 



UL 



lUl 



Figure 14. Wave patterns at consecutive time instants for parameters /ci = 1.0 and 
/c2 = l-l- The sohtons are moving in the same direction. 



6. Summary and discussion 

A simple classification scheme of the two soliton solutions of the Boussinesq equation 
have been presented. The scheme is based on the behavior of local maxima of the wave. 
We have shown that for solitons moving in the same direction there can be three different 
scenarios. In the Type I. case (see Figs. |3] and HI) there are two maxima all the time, 
separated by a minimum. In the Type Ila. case (see Figs. [5] and [6]) initially, when the 
solitons are still separated, there are two maxima with a minimum in between. During 
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the collision the two solitons merge and only one maximum remains, the remnant of the 
other shows up only as a shoulder. Later on, however, the second maximum reappears 
and grows. Then the first maximum disappears for a while, but as the solitons become 
separated, it reappears and we have the initial solitons in a reversed ordering along 
the line. In the Type lib. (see Figs. [7] and [H]) case the separated solitons merge to 
a wave having a single maximum, and later on a second maximum reappears and the 
solitons separate again. For solitons moving in opposite directions only the Type lib. 
case exists (see Figs. 19) and fTOl) . In that case the wave numbers ki and k2 may coincide 
(see Figs. [T5] and [T6|) . In contrast, for solitons moving in the same direction ki and k2 
must be different. If ki — )■ k2 (see Figs. [TSland [Ml) we have an extreme Type I. case, 
namely, there remains a large minimal distance between the solitons all the time, and 
the chasing soliton loses in height, while the escaping soliton gains in height during the 
collision. As a result, eventually the two solitons change their ordering along the line, 
without any close contact. 

Since our result are based on the long wave approximation, we expect that they 
should be observable in that limit. Also, the results for solitons moving in the same 
direction should follow from the Korteweg-de Vries equation as well, since that equation 
is obtained in the same approximation. 
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Figure 15. Intersections of the curve ([TT]) and ([T^ at identical parameters fci = ^2 
1.0 for solitons moving in opposite directions. 
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Figure 16. Wave patterns at consecutive time instants for identical parameters 
ki = k2 = 1.0. The solitons are moving in opposite directions. 
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